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ABSTRACT 

Supernovae (SNe) exploding in a dense circumstellar medium (CSM) are hypoth¬ 
esized to accelerate cosmic rays in collisionless shocks and emit GeV 7 rays and TeV 
neutrinos on a time scale of several months. We perform the first systematic search 
for 7 -ray emission in Fermi LAT data in the energy range from 100 MeV to 300 GeV 
from the ensemble of 147 SNe Type Iln exploding in dense CSM. We search for a 7 -ray 
excess at each SNe location in a one year time window. In order to enhance a possible 
weak signal, we simultaneously study the closest and optically brightest sources of our 
sample in a joint-likelihood analysis in three different time windows (lyear, 6 months 
and 3months). For the most promising source of the sample, SN 2010jl (PTFlOaaxf), 
we repeat the analysis with an extended time window lasting 4.5 years. We do not find 
a significant excess in 7 rays for any individual source nor for the combined sources 
and provide model-independent flux upper limits for both cases. In addition, we de¬ 
rive limits on the 7 -ray luminosity and the ratio of 7 -ray-to-optical luminosity ratio 
as a function of the index of the proton injection spectrum assuming a generic 7 -ray 
production model. Furthermore, we present detailed flux predictions based on multi¬ 
wavelength observations and the corresponding flux upper limit at 95% confidence 
level (CL) for the source SN2010jl (PTFlOaaxf). 

Keywords. Methods: data analysis; cosmic rays, gamma rays, supernova 


1. Introduction 


The Large Area Telescope (LAT) on-board the Fermi Gamma-ray Space Telescope mission 


unanticipatedly detected 7 -ray emission from five Galactic novae 


dAbdo et al.|20l0HHill et al.|2013 
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Hays et al. 2013; Cheung et al. 2013). The origin of the 7 -ray emission is still unclear. Shocks 


produced by expansion of the nova shell into the wind provided by the companion star or internal 
shocks within the ejecta might be responsible for acceleration of particles to relativistic energies and 
ensuing high-energy 7 -ray emission. A similar mechanism but with much larger energy output is 
hypothesized to produce 7 rays in supernovae (SNe) yielding potentially detectable 7 -ray emission 
even from extragalactic sources. Murase et al. ( |2011[ |2014[ ) and |Katz et ah] ( |2011| ) showed that 
if the SN progenitor is surrounded by an optically thick circumstellar medium (CSM), then a 
collisionless shock is necessarily formed after the shock breakout. The collisionless shock may 
accelerate protons and electrons to high energies, which emit photons from the radio-submillimeter 
through GeV energies and TeV neutrinos. Such conditions appear in shocks propagating through 
dense circumstellar matter (e.g., wind). Recently several candidates for such SNe powered by 
interactions with a dense CSM were found (e.g., Qfek et aL||2007> 2014a; Smith et al.||2009; Zhang 


et al. 

2012 

) and some superluminous supernovae were suggested to be powered by interactions 

(e.g., 

Quimby et al. 

2011 ; 

Chevalier & Irwin 

2011 

). Such interaction-powered supernovae may 


also be Pevatrons, implying their importance for the origin of the knee structure in the cosmic-ray 
spectrum (Murase et al. 2014; Sveshnikova 2003). Both 7 rays and neutrinos originate from pp 
and py interactions producing pions, which in the neutral case decay to 7 rays and in the charged 
case produce neutrinos in the decay chain. Thus, the initial neutrino and 7 -ray spectra have the 
same shape. Contrary to neutrinos, 7 rays might be affected by absorption in the CSM and/or 


two-photon annihilation with low-energy photons produced at the forward shock (Murase et al. 


2011). However, arguments made in Murase et al. (2014) suggest that GeV 7 rays can escape the 


system without severe attenuation, if the shock velocity is in the right range, especially late after 
the shock breakout. 


Motivated by the fact that the LAT has detected 7 -ray emission from novae, we are presenting 
the first systematic search for 7 -ray emission from Type Iln SNe in Fermi LAT data from 100 MeV 
to 300 GeV. Considering current theoretical uncertainties we are aiming for a model independent 
search. SNe positions and explosion times are given by optical surveys such as the Palomar 
Transient Factory (PTF - Law et al. 2009; Rau et al. 2009). 

We present the sample of SNe used in the 7 -ray data analysis in Section [ 2 } Section [3] de¬ 
scribes the Fermi LAT data analysis followed by an interpretation of our results in Section [4] and 
conclusions in Section [U 


Law et al. 2009; Rau et al. 2009 


2. SNe Sample 

Type Iln and Type Ibn SNe are the best candidates to be interacting with a dense CSM. Their 
long-lasting bright optical light curves are believed to be powered by the interaction of the ejecta 


with massive CSM (Svirski et al. 


2012 ). SNe of these types are often accompanied by precursor 
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mass-ejection events (Ofek et al. 2014b). Here we mainly use the PTF SN sample along with 


publicly available Type Iln SNe discovered since the launch of Fermi in 2008. Appendix [A] lists 
all 147 SNe of this sample that we consider in our 7 -ray search, i.e., all sources with an estimated 
explosion time later than 2008 August 4 and before 2012 May 1 (this is one year before the end of 
the studied 7 -ray data sample). The apparent A-band peak magnitude (to) as a function of the 
peak time is shown in Fig. [I] Note that throughout this paper we refer to m as the peak magnitude; 
for sources where the peak magnitude is not determined we use the discovery magnitude instead. 
The subsample of bright (to < 16.5) and/or nearby (with a redshift z < 0.015) SNe used for the 
joint likelihood analysis is detailed in Table [lj 



•yV° V tf .' 0 0 1' 0 7 ° &F 7 ° - S >'° 

date 


Fig. 1.— Apparent i?-band peak (detection) magnitude as a function of the peak (detection) time shown 
in red (blue) for all 147 SNe in our sample. For some SNe the peak time and magnitude is not determined; 
in those cases we use the detection time and magnitude. 
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Name 

RA (°)* 

Dec (°)* 

Date 

z 

m 

TS (p-value) 

SN2008gm 

348.55 

-2.78 

2008-10-22* 

0.012 

17.00* 

3.2 (0.169) 

SN2008ip 

194.46 

36.38 

2008-12-31* 

0.015 

15.70* 

0.0 (0.572) 

SN2009au 

194.94 

-29.60 

2009-03-11* 

0.009 

16.40* 

0.0 (0.572) 

PTFlOujc 

353.63 

22.35 

2009-08-05 

0.032 

16.20 

0.0 (0.572) 

SN2009kr 

78.01 

-15.70 

2009-11-06* 

0.006 

16.00* 

4.7 (0.104) 

SN2010bt 

192.08 

-34.95 

2010-04-17* 

0.016 

15.80* 

14.4 (0.0065) 

PTFlOaaxf 

SN2010jl 

145.72 

9.50 

2010-11-18 

0.011 

13.20 

7.1 (0.039) 

PTFlOaaxi 

SN2010jp 

94.13 

-21.41 

2010-11-23 

0.010 

18.00 

0.0 (0.572) 

SN2011A 

195.25 

-14.53 

2011-01-02* 

0.009 

16.90* 

0.0 (0.572) 

PTFlliqb 

8.52 

-9.70 

2011-08-06 

0.013 

15.20 

0.3 (0.469) 

SN2011fh 

194.06 

-29.50 

2011-08-24* 

0.008 

14.50* 

1.9 (0.262) 

PSNJ10081059+5150570 

SN2011ht 

152.04 

51.85 

2011-10-29 

0.004 

14.50 

0.0 (0.572) 

PTFllqnf 

86.23 

69.15 

2011-11-01* 

0.014 

19.80* 

1.4 (0.320) 

SN2011hw 

336.56 

34.22 

2011-11-18’!' 

0.023 

15.70* 

0.0 (0.572) 

SN2012ab 

185.70 

5.61 

2012-01-31* 

0.018 

15.80* 

0.0 (0.572) 

PSNJ18410706-4147374 

SN2012ca 

280.28 

-41.79 

t 

2012-04-25 

0.019 

t 

14.50 

0.0 (0.572) 


* Epoch J2000.0 

* Discovery date 

* Discovery magnitude 


Table 1: List of nearby and/or bright SNe - with redshift z < 0.015 and/or 11- 
band magnitude m < 16.5. The colums contain the name of the SN, its direction 
in equatorial coordinated (right ascension, RA, and declination, Dec), its peak 
date and peak R-band magnitude, its redshift, its test statistic (TS) and p-value. 
See Section B for details on the TS and p-value calculation. Note that if the 
peak date and magnitude are not available in the catalog, the discovery date and 
magnitude are quoted instead. 
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3. Fermi LAT 7 -ray Data Analysis 


The Fermi LAT is a pair-conversion telescope, sensitive to 7 rays with energies from 20 MeV 


to greater than 300 GeV (Atwood et al. 2009). It has a large field of view and has been scanning 


the entire sky every few hours for the last 6 years. Thus it is very well suited for searches 
for transient 7 -ray signals on the timescale of months. Bright SN events may be detectable at 


distances d < 30Mpc (Murase et al. 2011) depending on the properties of the source. Margutti 


et al. (2014) searched for 7 rays from a single SN in the case of SN2009ip located at a distance of 


24Mpc. No 7 -ray excess was identified in Fermi LAT data at the SN position; this is consistent 
with the picture of ejecta colliding with a compact and dense but low-mass shell of material. 
For a detection of a single source, closer and/or brighter SNe are needed (i.e., reached by larger 
dissipation and larger CSM masses). The properties of the ejecta and CSM can be estimated from 
multi-wavelength observations in a few cases (e.g., SN2009ip, Ofek et al. 2013), but are uncertain 
or not known in most cases. 

In this analysis we use 57 months of Fermi LAT data recorded between 2008 August 4 and 2013 
May 1 (Fermi Mission Elapsed Time 239557418-389092331 s), restricted to the Pass 7 Reprocessed 
Source clasf[/j We select the standard good time intervals (e.g., excluding time intervals when the 
field of view of the LAT intersected the Earth). The Pass 7 Reprocessed data benefit from an 
updated calibration that improves the energy measurement and event-direction reconstruction 
accuracy at energies above 1 GeV (Bregeon et al. 2013). To minimize the contamination from the 
7 rays produced in the upper atmosphere, we select events with zenith angles < 100°. We perform 
a binned analysis (i.e., binned in space and energy) using the standard Fermi LAT ScienceTools 
package version v09r32p05 available from the Fermi Science Support Centeij^] (FSSC) using the 
P7REP_SOURCE_V15 instrument response functions. We analyze data in the energy range of 
100 MeV to 300 GeV binned into 20 logarithmic energy intervals. For each source we select a 
20° x 20° region of interest (ROI) centered on the source localization binned in 0?2 size pixels. 
The binning is applied in celestial coordinates and an Aitoff projection was used. 


We use four different approaches in our analysis: 


1. We perform a likelihood analysis to search for 7 -ray excesses that are consistent with origi¬ 
nating from a point source coincident with the position of each Type Iln SNe in our sample 
over a 1 -year time scale. We assume that their 7 -ray emission follows a power-law spectrum. 
This approach is sensitive to single bright sources. 

2 . In a model-independent approach (i.e., no prior assumption on the SN 7 -ray spectral shape) 
we compute the likelihood in bins of energy (bin-by-bin likelihood). We use the bin-by-bin 


J http://fermi.gsfc.nasa.gov/ssc/data/analysis/documentation/Pass7REP_usage.html 
2 http://fermi.gsf c.nasa.gov/ssc/data/analysis/ 
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likelihood to evaluate 95% CL flux upper limits in 20 energy bins for the 16 closest and 
optically brightest SNe in our sample. 

3. In order to increase the sensitivity for a weak signal, we combine individual sources in a joint 
likelihood analysis using the composite likelihood tool, Composite2 , of the Fermi Science 
Tools. 

4. We repeat the joint likelihood analysis using the composite likelihood tool, but limit the 
sample to those Type Iln SNe that exhibit additional indications of strong interactions with 
their CSM. Not all Type Iln SNe might be surrounded by a massive CSM. This clean sample 
of SNe with confirmed massive CSM might produce a strong 7 -ray signal and should provide 
an enhanced signal-to-background ratio. 


Accurate SN positions are given by optical localizations. Theoretical predictions of the duration 
of the 7 -ray emission are uncertain and motivate a search in several time windows. We test 
three different time windows: AT = lyear, 6 months and 3 months. The optical light curve is 
produced by the interaction of the SN ejecta with the dense CSM and is thus correlated with 
the expected 7 -ray emission. Most of the 7 -ray emission is expected during the interactions after 
the shock breakout. The optical light curve peak is reached around the end of the breakout (see 
e.g., Ofek et al. 2010 ). We collected the SN properties from the PTF sample, Astronomer’s 
Telegrams (ATels^and Central Bureau for Astronomical Telegrams (CBETsj^} Most PTF sources 
are unpublished and the other events were drawn from ATEL and CBET. Full details and final 
analysis of the PTF SN Iln sample will be provided in a forthcoming publication. In some cases 
the known SN properties include the optical flux peak time while in other cases this information 
is missing and only the optical detection time is available. To account for the uncertainty in the 
determination of the peak time and to make sure no early 7 -ray emission is missed, we start the 
time window 30 days before the peak time (or the detection time in case the peak time is not 


provided). In the case of the three novae, the reported 7 -ray light curves (see Fig. 1 in Hill 


[et al. 2013) have very similar durations justifying a similar time window for all sources. However, 
the duration of the novae detected by Fermi were ~ 20 days, while SNIIn typically last longer, 
(9(100days - lyear). 


3.1. Source Specific Analysis 

We analyze the 20° x 20° ROI around each source in our SN sample in a 1 -year time window 
in a binned likelihood analysis. We construct a model whose free parameters are fitted to the 

3 http://www.astronomerstelegram.org/ 

4 http://www.cbat.eps.harvard.edu/cbet/RecentCBETs.html 
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data in the ROI. This model includes a point-like source at the SN position; its 7 -ray spectrum is 
represented as a power-law function with both index and normalization free to vary. In addition 
we have to model the point sources in the ROI and the diffuse 7 -ray emission. We consider all 
the 2 FGL sources (Nolan et al. 2012 ) included within a larger region of radius, R = 20°, to 
allow for the breadth of the LAT point-spread function that may cause a significant signal from 
sources outside the ROI to leak into it. The positions and spectral parameters of all 2 FGL sources 
within 15° < R < 20° from the center of the ROI are fixed to the values reported in the 2 FGL 
catalog; those are on average 21 sources. For the sources within 5° < R < 15° with > 15<r 
detection significance in 2 FGL only the flux normalization is let free to vary and all the other 
parameters are fixed to the values reported in the 2 FGL catalog. The parameters for all the other 
sources within 5° < R < 15° are fixed to the 2 FGL catalog values. Finally, for sources within 
R < 5° all parameters (index and normalization in case of a power-law spectrum, index, cutoff 
and normalization in case of a power-law with exponential cutoff and normalization, spectral slope 
and curvature in case of a log-parabola source spectrum) are free to vary if the source significance 
exceeds 4 a, otherwise all source parameters are fixed. On average 3 sources per ROI have all 
parameters free, while 6 sources have a free normalization and 18 sources are fixed to the 2 FGL 
values. 


We determine the best values for all the free parameters fitting our source model together 
with a template for the isotropic and Galactic interstellar emission)^ to the LAT data with a binned 
likelihood approach as described in Abdo et al. (2009). To quantify the significance of a potential 
excess above the background, we employ the likelihood-ratio test (Neyman & Pearson 1928). We 
form a test statistic 

TS = — 2 Alog£ = — 2 (log £ 0 — log£), ( 1 ) 


where £ 0 is the likelihood evaluated at the best-fit parameters under a background-only, null 
hypothesis, i.e. a model that does not include a point source at the SN position, and C the 
likelihood evaluated at the best-fit model parameters when including a candidate point source at 
the SN position. 


The distribution of the TS values obtained for all the SNs using a 1 year time window is 
displayed in Fig. [ 2 ] (left), compared to the TS distribution obtained from performing a similar 
analysis at random positions in the sky. We require the random ROI centers to be separated by 
at least 3.5° and to lie outside of the Galactic plane, i.e. | 6 | > 10°. The analysis in the Galactic 
plane region is complicated by the intense Galactic diffuse emission and none of the SNe in our 
sample is located close to the plane. Those requirements limit the number of independent ROIs; 
we use 1140 ROIs in our analysis. The distribution of SN-position TS values is similar to the 


5 We use the templates provided by the FSSC for the P7REP_SOURCE_V15 event class (http://fermi.gsfc. 
nasa.gov/ssc/data/access/lat/BackgroundModels.html) with free normalization and free index in case of the 
Galactic interstellar emission model. 
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distribution of random-position TS values (see Fig. [ 2 ] left). The highest TS value found among 
the SN positions is 14.4, which corresponds to a p-value of 0.0065 (obtained from the random 
position analysis), which is below 3 a for a single trial (see Fig. fright). Given the number of SNe 
in our sample a trials factor needs to be applied, which increases the p-value to 0 . 6 . 

Optically bright SNe are expected to produce a brighter 7 -ray signal than optically dim ones 
and nearby SNe are expected to be brighter than sources at large distance. However, we do not 
find an obvious correlation of TS value with redshift or magnitude (see left and right panels of 
Fig-! respectively), indicating that the 7 -ray signals of individual SNe, if present, are weak. 

Three of 147 SNe have a 2 FGL source in their close vicinity with an angular distance of less 
than 0.4°. In each case the nearest 2 FGL source is associated with an active galactic nucleus 
through multi-wavelength data. Since the spectral parameters of the nearby source are left free 
to vary in the fit, a possible SNe flux could have been absorbed by the background source. Those 
sources are PTFlOweh, LSQ12by and SN 2012bq, which are optically dim and distant sources and 
thus not part of the subsample of nearby and/or bright SNe. 




Fig. 2.— Left: Distribution of TS values for a test source modeled by a power-law energy spectrum 
located at a SN position (red), compared to TS for a similar test source located at a random position 
(grey). Right: Cumulative distribution of random-position TS values. The blue dashed (dotted) line 
indicates a Gaussian equivalent one-sided 3 a (2 a) probability of finding a larger TS than the TS indicated 
by the intersection of the blue line with the grey distribution. The red solid line shows the largest TS 
found in the source-specific analysis, which has a p-value of 0.0065 (red dotted line) and thus lies below 
3 a. Considering the trials factor, the p-value increases to 0.6. 


3.2. Model independent analysis of nearby and/or bright SNe 

The 7 -ray spectral shape resulting from particle acceleration in the interaction of SN ejecta 
with a dense CSM is not known a priori. It is determined by the initial proton spectrum and 
could be altered by absorption of the 7 rays in the surrounding medium. Therefore, we study 
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Fig. 3.— TS of a test source located at the SN position as function of redshift (left) and as function of 
magnitude (right). Note that the largest TS value was 14.4 (corresponding to a p-value of 0.009), which 
does not exceed the detection threshold of 5 a. 


the closest and/or optically brightest sources, which are the most promising sources in terms of 
expected 7 -ray emission, in an approach independent of a SN spectral model assumption. The 
sources chosen for this analysis have to fulfill the criteria of 2 < 0.015 or m < 16.5, and are listed 
in Table [l] We fix the spectral parameters of the background sources and the diffuse templates 
to their global values obtained from the source-by-source analysis over the entire energy range 


described in Section 3.1 (without including the SN itself). Following the procedure described 
in Ackermann et al. (2014) we calculate the likelihood in each of the 20 energy bins after inserting 


a test source at the SN position at various flux normalization values: 


( 2 ) 


where Vj is the photon data, C :I the Poisson likelihood and {/u,} a set of independent signal 
parameters in energy bin j. The symbol 0 represents the nuisance parameter (i.e., free parameters 
of background sources and diffuse templates) and 6 indicates that they have been fixed to their 
global values. The bin-by-bin likelihood allows us to find the upper limits at 95% Clfl defined as 
the value of the energy flux, where the log-likelihood decreases by 2.71/2 from its maximum (the 
“delta-log-likelihood technique” - [Bartlett 1953; Rolke et al. 2005). An example is shown in 
Fig. [i] for SN 2010J1 (PTFlOaaxf - |Ofek et al. 2014a; |Fransson et al.||2014 ; Zhang et al. 2012), 
while similar plots for all nearby sources can be found in Appendix [B] Any SN model predicting 
a certain 7 -ray spectrum can be tested using those results (see Ackermann et al. 2014, for more 
details on the bin-by-bin likelihood) by recreating a global likelihood by tying together the signal 
parameters over the energy bins: 


£ 0,0 I'D) = 


( 3 ) 


6 Note, we are using a two-sided confidence interval. 
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with n denoting the global signal parameters. 


For the most promising source of our sample, SN 2010jl, we repeat the analysis for an extended 
time window ending in May 2015, i.e., spanning 4.5 years. This is motivated by the fact that in 


some cases SN Type Iln emission lasts for 3-5 years after the explosion (Cooke et al. 2009) 
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Fig. 4.— Histogram of the bin-by-bin LAT likelihood function used to test for a putative 7 -ray source 
at the position of supernova SN 2010jl (PTFlOaaxf). The bin-by-bin likelihood is calculated by scanning 
the integrated energy flux of the SN within each energy bin (equivalent to scanning in the spectral 
normalization of the source). When performing this scan, the flux normalizations of the background 
sources are fixed to their optimal values as derived from a maximum likelihood fit over the full energy 
range. Within each bin, the color scale denotes the variation of the logarithm of the likelihood with 
respect to the best-fit value of the SN flux using a 1 year time window. Upper limits on the integrated 
energy flux are set at 95% CL within each bin using the delta-log-likelihood technique and are largely 
independent of the SN spectrum. The black arrows indicate the 95% CL flux upper limits for AT = 1 yr, 
where the shown log-likelihood decreases by 2.71/2 from its maximum. For completeness we overlay the 
95% CL upper limits for AT = 6 months and AT = 3 represented by dotted-dashed and dotted lines 
respectively. For the particular case of SN 2010jl we repeated the analysis for an extended time window 
spanning 4.5 years. The results are overlaid as cyan dotted line. 


3.3. Joint Likelihood Analysis 


For greater sensitivity to a weak 7 -ray signal from interaction-powered SNe, we combine the 
16 closest and/or brightest sources in a joint likelihood analysis. To be independent from any 


spectral shape assumption we perform the analysis in energy bins (see Section 3.2 for details of 
the bin-by-bin likelihood analysis). In each energy bin we tie the SN flux normalization for all 
16 SNe together resulting in one free parameter per energy bin. The likelihood values for the 
individual sources, i, are multiplied to form the joint likelihood: 


c^mv) = \[c^Am 


( 4 ) 
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However, we have to make some assumption about a common scaling factor of the 7 -ray flux 
in order to tie the SNe flux normalizations together (i.e., we want to give a larger weight to 
SNe with greater expected 7 -ray fluxes in the joint likelihood). We use two different approaches: 
first, we assume that all SNe have the same intrinsic 7 -ray luminosity; therefore, the observed 
7 -ray flux for each SN scales with a factor inversely proportional to the square of the luminosity- 
distance d. The redshift is measured for each SNe and since we only consider nearby SNe we 
use a simple linear approximation for the relation between redshift and distance: d = z x c/H , 
with H = 67.8kms _1 Mpc -1 (Ade et al. 2014). We do not apply a redshift-dependent energy 


rescaling for SNe at different redshifts, since the energy shift is negligible at the small redshifts 
(i.e., z < 0.015) considered in this analysis. We weight the flux normalization in each energy bin 
of each source with uq = (lOMpc/d) 2 . We then tie those weighted normalizations together. The 
exact value of H does not influence our results since the combined normalization of all sources is 
free in the fit of the model to the data in each energy bin. Note that only the SN flux normalization 
is free while the background source parameters as well as the diffuse template parameters are fixed 
to their global values obtained from a fit to the entire energy range. 

Alternatively, we assume that the 7 -ray flux is correlated with the optical flux, i.e. we use a 
weight proportional to the optical flux^Jor 10 -a4m . We chose the weight to be: 


w ri 


_ 2Q—0.4(m— C) 0-4tti+5.2 


( 5 ) 


where m is the apparent A-band magnitude provided by the SN catalog and C = 13 is a normal¬ 
ization constant. Again, the exact choice of C does not influence our results since the combined 
normalization of all sources is free in the fit. We chose to neglect a correction for Galactic dust 
extinction, which is at most 0.28 mag and thus smaller than the uncertainty in the peak magnitude 
determination. 

We perform the joint likelihood analysis for three time windows: lyear, 6 months, 3 months 
since the A-band maximum light. Fig. [5] shows the likelihood profiles of the combined 7 -ray flux. 
Table [ 2 ] summarizes the results from the combined likelihood analysis and shows the sum of TS 
over all energy bins. No significant improvement in the likelihood by including the SNe in the fit 
could be found in the joint likelihood analysis. The largest TS value of 8.8 is found in case of 
assuming the 7 -ray flux scales with the optical flux for the 1 -year time windows. According to 
Wilks’ theorem, TS is distributed approximately as x 2 with the degrees of freedom equal to the 
number of parameters characterizing the additional source. Taking into account the number of 
free parameters ( 20 , one for each energy bin) the probability that this is a statistical fluctuation 
is 98.5%. This significance would be further decreased by taking into account trials factors for the 
two different weighting schemes and 3 different time windows. 


7 Note that flux and apparent magnitude are related through: m — mo = — 2.51og 10 where Fo and too are 

the flux and apparent magnitude of a reference star. 
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However, if we assume a spectral model for the SN flux, we can greatly reduce the number 
of free parameters. For illustration we fit a power-law spectral shape to the bin-by-bin likelihood 
following Equation [3| The index and normalization of the power-law function are left free to vary 
in the fit. The resulting TS values and corresponding p-values (not including trials factors) are 
summarized in Table [2j none of them are significant. A more physical spectral model is fitted to 
the bin-by-bin likelihood in Section |4} 



Fig. 5.— Similar to Fig. [4| but for the composite likelihood instead of the single-source likelihood. 
Left: Composite likelihood profile for each energy bin weighting each source with (10Mpc/d) 2 . Right: 
Composite likelihood profile for each energy bin weighting each source with io -a4,n + 5 - 2 . The black arrows 
indicate the 95% upper limits for AT = 1 yr, while the dotted-dashed and dotted lines represent the 95% 
upper limits for AT = 6 months and AT = 3 months, respectively. 


Weighting TS TS pi (p-value) 



1 yr 

6 months 

3 months 

lyr 

6 months 

3 months 

(10Mpc/d) 2 

2.2 

2.1 

2.4 

0.0 (1.0) 

0.0 (1.0) 

0.0 (1.0) 

IQ— 0.4ra+5.2 

11.7 

7.8 

9.0 

2.9 (0.23) 

1.6 (0.45) 

0.0 (1.0) 


Table 2: Sum over bin-by-bin TS values obtained from the joint likelihood analysis. TSpp is the 
TS obtained by assuming a power-law spectral shape. 


3.4. Joint Likelihood Analysis of SN Subsample with Confirmed Massive CSM 

We select a subsample of 16 SNe from the Type Iln SNe catalog for which we have additional 
evidence through multi-wavelength observations for the existence of a massive CSM. We select 
SNe that show Balmer emission lines and continuum in both early and late times. The SNe in 
this sample are: PTF12csy, PTFlloxu, PTFllmhr, PTFllfzz, PTFllfuu, PTFlOaaxf, PTFlOptz, 
PTFlOscc, PTFlOjop, PTFlOfei, PTFlOqaf, PTFlOtel, PTFlOtyd, PTFlOgvf, PTFlOcwl, PTF09drs 
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We repeat the joint likelihood analysis described above for this subset with the optical flux weight¬ 
ing scheme for three time windows (1 year, 6 months and 3 months). The results are displayed in 
Fig-i The TS values of the composite fit are 11.3, 17.5 and 10.3 for the time windows of 1 year, 
6 months, and 3 months, respectively. Taking into account the 20 free parameters, the chance 
probability for a TS of 17.5 is 62%. 
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Fig. 6.— Joint likelihood analysis of the SN subsample with confirmed massive CSM: Joint likelihood 
profile for each energy bin weighting each source with 10~°' 4m+5 ' 2 . The black arrows indicate the 95% 
upper limits for AT = 1 year, while the dotted-dashed and dotted lines represent the 95% upper limits 
for AT = 6 months and AT = 3 months, respectively. 


4. Interpretation 


Murase et al. (2011) suggested that 7 -ray emission is produced by cosmic rays accelerated 


at the early collisionless shock between SN ejecta and circumstellar material. For the scenario 


described by Murase et al. (2014), 7 -ray emission can be predicted when the model parameters are 


determined by optical and X-ray observations. We defer such model-dependent analyses to future 
work. Instead, in this work, we take a model-independent approach, where we aim to constrain the 
7 -ray luminosity as a function of the proton spectral index. We assume that the spectrum of CR 
protons is given by a power law (in momentum) with the minimum and maximum proton momenta 


of 0.1 GeV/c and 10 8 GeV/c, respectively. Then, we calculate the 7 -ray flux following Kelner et al. 
(2006). In the calorimetric limit, which is expected for SNe like SN 2010jl (Murase et al. 2014), 
the 7 -ray spectral index follows the proton spectral index, although the resulting limits (shown in 
Fig.§ are similar to what would be obtained for non-calorimetric cases, for which the resulting 
shape of the 7 -ray spectrum is slightly harder than the proton spectral shape due to the energy 
dependence of the pp cross section. For simplicity we do not take into account 7 -ray absorption; 


Murase et al. (2014) showed that GeV 7 rays can escape from the system without severe matter 


attenuation if the shock velocity is high enough. 
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The diffusive shock acceleration theory predicts that the proton acceleration efficiency is 
0.1. In the calorimetric limit, all the proton energy is used for pion production, and 1/3 of 
pions are neutral pions that decay into 7 rays. Then, about half of the 7 rays are absorbed deep 
inside the ejecta, so we expect L 1 « (l/ 6 )e ; ,/ csc Li dll , where L\^ u is the kinetic luminosity and / esc 
is the escape fraction of 7 rays. The 7 -ray attenuation due to the Bethe-Heitler process is relevant 
when the shock velocity is lower than ~ 4500 km s _1 , while the two-photon annihilation process 


6 p ~ 


is relevant when the shock velocity is high enough (Murase et al. 2014). Although 7 rays can 


escape late after the shock breakout, the attenuation can be relevant around the shock breakout 
so we assume f esc ~ 0.1-1 to take into account uncertainty of the 7 -ray flux. The radiation energy 
fraction is given by e 7 = L ra( j / Li dn . where L Ta( j is the bolometric radiation luminosity. About 
half of the kinetic energy is converted into the thermal energy, and half of the thermal energy is 


released as outgoing radiation, which implies e 7 ~ 1/4 (Ofek et al. 2014a). As a result, we have 
L 7 /L ra d ~ (l/6)(e p /e 7 )/esc ~ (l/15)/ esc . Our limits presented below are on the fraction of 7 -ray 
to R-band luminosity, which is an upper bound on L 7 /L ra d. In the case of SN 2010J1 Lr ~ L ra< i 
and thus L 1 /Lr ~ 0.01-0.1 is theoretically expected. 

As an example, we consider supernova SN 2010J1 (PTFlOaaxf), which is the most-likely de¬ 
tectable CR accelerator, because multi-wavelength observations indicate a very massive CSM 
of 10 M q . We present a generic flux prediction for the calorimetric limit for this source as¬ 
suming a proton spectral index of T p = —2 and a normalization of the 7 -ray flux that yields 
0.01 < L 1 /Lr <0.1 (shown as shaded green region in Fig. [ 7 ]) and calculate the corresponding 
flux upper limit (shown in blue in Fig. [7|) following the procedure outlined in Ackermann et al. 


(2014). The bin-by-bin likelihood analysis is used to re-create a global likelihood for a given signal 


spectrum by tying the signal parameters across the energy bins (see Eqn. [3j). In this case the 
global signal parameter is the flux scale factor N relative to the flux that yields L~,/Lr = 0.1 (i.e., 
the upper bound of the uncertainty band shown in Fig. [Tj left). We assume that SN 2010jl is at 
distance 48.7 Mpc with an apparent R-band peak magnitude of 13.2. We calculate the change in 
log-likelihood for various values of N and find the 95% flux upper limit (given by the value of N for 
which the delta log-likelihood decreases by 2.71/2 compared to its minimum). The derived upper 
limit touches the optimistic model prediction, i.e. the upper bound of the theoretical uncertainty 
band. A more detailed modeling of the expected flux based on multi-wavelength observations is 
outside the scope of this paper and will follow in future work. Better constraints on the 7 -ray 
escape fraction are crucial to calculate stringent limits on the proton acceleration efficiency and 
will be obtained in more detailed modeling. 

More stringent limits are expected from the joint likelihood result^ Generic 7 -ray flux 
predictions for various proton spectral indices are shown in Fig. [ 8 } We calculate the 95% CL 


8 Note that including sources with a statistical over-fluctuation can worsen the joint limit. 
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Fig. 7.— Left: Predicted 7 -ray energy spectrum for SN 2010jl assuming T p = — 2 and a normalization of 
the 7 -ray flux yielding 0.01 < L 1 /Lr <0.1 shown as the green shaded region compared to the 95% flux 
upper limit (blue). Right: Likelihood profile for the spectral normalization parameter N relative to the 
flux prediction yielding L 1 /Lr = 0.1. The dashed green line indicates an increase of the negative delta 
log-likelihood by 2.71/2 compared to its minimum. 


upper limit on the 7 -ray luminosity 

F 1 

— AttcPF^ = 47 r (10 Mpc ) 2 —, 
7 w d 


( 6 ) 


where F/ is the integrated 7 -ray flux over the energy range used in this analysis. The luminosity 
L~, is proportional to the result of the joint likelihood analysis using the weight wa = (lOMpc/d) 2 , 
assuming all sources have the same L 7 . In other words our joint likelihood results set a limit on 
F^/wd and thus on L 1 . The result is shown in Fig. [ 9 ] (left) as a function of the proton spectral 
index. 


In addition we calculate the 95% CL upper limit on the ratio of 7 -ray to optical luminosity 

47 rd 2 < _ 47t(1 Mpc ) 2 F* 

7/ R L o 10°- 4 ( M ®- M ) Z/ O 10 °- 4M © +4 - 8 w m ’ 1 1 

where L 0 = 6 x 10 32 ergs -1 is the R-band luminosity and M 0 = 4.7 the absolute R-band magnitude 
of the Sun. The ratio is proportional to F 1 /w m , which is constrained by the joint likelihood analysis 
assuming a correlation of optical and 7 -ray flux, i.e. weighting with w m = 10 -a4m+52 . Thus we 
can use the joint likelihood results to set a limit on L 1 /Lr as a function of (see Fig. fright). 

In Fig. [9] both limits discussed above are compared to the limit obtained using only one SN. 
The closest SN (SN2011ht with a distance of d — 17.7Mpc) in the case of 1/d 2 weighting and 
the brightest SN (SN 2010jl with a magnitude of m = 13.2) in case of weighting with the optical 
flux. In both cases the combined limit is dominated by one SN. In the case of 1/d 2 weighting the 
single source limit is better than the combined limit, indicating a statistical under-fluctuation in 
the individual analysis of this source or an over-fluctuation in one of the sources included in the 
joint likelihood. 
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Fig. 8 .— Gamma-ray energy spectra assuming a total 7 -ray energy of 10 50 erg, a source distance of 
lOMpc and a duration of 1 year for various proton spectral indices T p . The shaded grey region shows 
the energy range covered by this analysis. 




Fig. 9.— Left: 95% CL upper limit on the 7 -ray luminosity as a function of the proton spectral 
index based on the results obtained from the joint likelihood analysis with 1/d 2 weighting shown in blue 
compared to the limit obtained from the closest single source SN 201 lht in green. Right: 95% CL upper 
limit on the ratio of 7 -ray and optical luminosity L^/Lr as a function of T p assuming a proportionality 
between optical and 7 -ray flux shown in blue compared to the limit obtained from a single source analysis 
of SN 2010jl considering a one year time window (in green). The results of the analysis with an extended 
time window of 4.5 years for SN 2010jl are shown in dashed green. 
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5. Conclusions 

The origin of the multi-wavelength emission of Type Iln SNe and the onset of cosmic-ray 
production in supernova remnants is not fully understood. Type Iln SNe are expected to be host 
sites of particle acceleration, which could be pinpointed by transient 7-ray signals. For the first 
time we searched in a systematic way for 7-ray emission from a large ensemble of Type Iln SNe in 
coincidence with optical signals. No evidence for a signal was found, but our observational limits 
start to reach interesting parameter ranges expected by the theory. We set stringent limits on 
the 7-ray luminosity and the ratio of 7-ray and optical luminosity. For example, we can exclude 
L 7 /Lr > 0.1 at 95% CL for proton spectral indices of < 2.7 from the results of the combined 
likelihood analysis assuming that L 1 /L R is constant. Those constraints can be converted to limits 
on the proton acceleration efficiency. In the case of SN 2010jl, our limits are close to theoretically 
expected values. However, uncertainties in the modeling, including the 7-ray escape fraction, leads 
to the range of 0(10%) to 0(1%) for the ratio of 7-ray to optical luminosity. Model-dependent 
calculations based on multi-wavelength observations will be performed in a future work and will 
allow us to set stringent constraints on the proton acceleration efficiency. 

We do not have to make this assumption in the analysis of individual SNe. The results from 
the optically brightest SN in our sample, SN 2010jl, alone leads to only a factor of two weaker 
constraints, excluding L 7 / L R > 0.2. Assuming a scaling of the 7-ray flux with 1 / d 2 we can exclude 
L 1 > 4 x 10 40 erg s _1 at 95% CL for all indices considered. A total 7-ray luminosity of 10 50 erg 
emitted within 1 year (as assumed in Fig. [8j) is excluded. The limits presented here are based on 
minimal assumptions about the 7-ray production and can be used to test various models. 
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A. SN Catalog 

The following table contains all SNe included in this analysis. The column definition is similar 
to Table [H 


Name 

RA (°)* 

Dec (°)* 

Date 

z 

m 

TS (p-value) 

SN2008gm 

348.55 

-2.78 

2008-10-22* 

0.012 

17.00* 

3.2 (0.169) 

CSS081201-103354-032125 

158.47 

-3.36 

2008-12-01* 

0.060 

18.30* 

0.0 (0.572) 

CSS080701.234413+075224 

SN2008ja 

356.05 

7.87 

2008-12-30 

0.069 

18.50 

0.0 (0.572) 

SN2008ip 

194.46 

36.38 

2008-12-31* 

0.015 

15.70* 

0.0 (0.572) 

SN2009au 

194.94 

-29.60 

2009-03-11* 

0.009 

16.40* 

0.0 (0.572) 

CSS080928_160837+041626 

SN2008iy 

242.16 

4.27 

2009-03-21 

0.041 

17.60 

0.3 (0.458) 

SN2009cw 

226.26 

48.67 

2009-03-28* 

0.150 

20.30* 

0.0 (0.572) 

SN2009eo 

224.53 

2.43 

2009-05-13* 

0.044 

18.10* 

0.0 (0.572) 

SN2009fs 

274.80 

42.81 

2009-06-01* 

0.054 

17.00* 

3.5 (0.154) 

PTF09ij 

218.06 

54.86 

2009-06-03 

0.124 

20.30 

0.0 (0.572) 

PTF09ge 

224.26 

49.61 

2009-06-06 

0.064 

17.90 

3.3 (0.165) 

PTF09tm 

206.73 

61.55 

2009-06-25 

0.034 

16.80 

0.0 (0.572) 

PTF09uj 

215.05 

53.56 

2009-06-26 

0.066 

18.20 

0.0 (0.572) 

PTF09uy 

190.98 

74.69 

2009-07-03 

0.313 

19.40 

0.0 (0.572) 

PTF09bcl 

271.61 

17.86 

2009-07-19* 

0.062 

20.87* 

0.0 (0.572) 

PTFlOujc 

353.63 

22.35 

2009-08-05 

0.032 

16.20 

0.0 (0.572) 

PTF09drs 

226.63 

60.59 

2009-08-15 

0.045 

18.50 

0.0 (0.561) 
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CSS090925.001259+144121 

3.25 

14.69 

2009-09-25+ 

0.090 

18.80+ 

0.0 (0.568) 

SN2009ma 

127.24 

0.59 

2009-10-17+ 

0.089 

18.20+ 

0.0 (0.572) 

CSS091018-091109+195945 

SN2009mb 

137.79 

20.00 

t 

2009-10-18 

0.150 

t 

19.00 

0.0 (0.572) 

SN2009kn 

122.43 

-17.75 

2009-10-26+ 

0.016 

16.60+ 

0.0 (0.572) 

SN2009kr 

78.01 

-15.70 

2009-11-06+ 

0.006 

16.00+ 

4.7 (0.104) 

SN2009nm 

151.35 

51.28 

2009-11-20+ 

0.210 

18.80+ 

0.0 (0.572) 

CSS091217-110637+341952 

SN2009nj 

166.65 

34.33 

t 

2009-12-17 

? 

1 

18.70 

0.0 (0.572) 

CSS091218-104011+223735 

SN2009nw 

160.05 

22.63 

t 

2009-12-18 

0.140 

t 

19.40 

0.0 (0.572) 

PTFlOdk 

77.09 

0.21 

2009-12-18+ 

0.074 

20.14+ 

0.0 (0.572) 

PTFlOu 

152.49 

46.01 

2010-01-05 

0.150 

19.80 

0.0 (0.572) 

PTFllner 

125.58 

72.83 

2010-01-11+ 

0.117 

20.94+ 

0.0 (0.572) 

PTFlOts 

SN2009nn 

188.49 

13.92 

2010-01-12 

0.046 

17.66 

7.9 (0.033) 

CSSlOOl 13.032138+263650 

SN2010M 

50.41 

26.61 

t 

2010-01-13 

0.060 

t 

18.80 

0.1 (0.517) 

PTFlOcwl 

189.09 

7.79 

2010-03-13+ 

0.085 

19.00+ 

0.0 (0.572) 

SN2010al 

123.57 

18.44 

2010-03-13+ 

0.017 

17.80+ 

9.7 (0.023) 

PTFlOcwx 

188.32 

-0.05 
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B. Likelihood Profiles in Energy Bins 


Figures |10| and |11| show the likelihood 
upper limit for three time windows AT = 
listed in Table [lj 


profiles in energy bins for 
1 yr, AT = 6 months and 


AT 

AT 


1 yr and the 95% CL 
3 months for all SNe 


9 This source is of Type Ibn, while all other sources are of Type Iln. 
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Fig. 10.— Similar to Fig. [4| Colors represent the likelihood profile for each energy for AT = lyr. The 
black arrows indicate the 95% CL upper limits for AT = 1 yr, while the dotted-dashed and dotted lines 
represent the 95% CL upper limit for AT = 6 months and AT = 3 months respectively. 


—Alog£ — Alog£ — Alog£ -Alog£ 
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Fig. 11.— Similar to Fig. [4| Colors represent the likelihood profile for each energy for AT = lyr. The 
black arrows indicate the 95% CL upper limits for AT = 1 yr, while the dotted-dashed and dotted lines 
represent the 95% CL upper limit for AT = 6 months and AT = 3 months respectively. 


—Alog£ — Alog£ — Alog£ -Alog£ 
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